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16  June,  discussions  were  held  at  MIRADCOM  with  Dwight  McPherson 
(Boeing)  ahd  Robert  Smith  (NWC/Chlna  Lake).  As  a  result  of  these 
discussions,  two  revisions  have  been  made  to  the  extended  target  simulation 
program!^ 

target  aspect  angle  Is  measured  from  the  roll  axis, 
not  In  the  wing  plane.  This  change  was  Implemented  In 
ETGEO  by  redefining  a  new  angle  (ANGL),  and  In  the  argument 
_ list  In  SCTAMP^  ^  (P 

In  order  to  Introduce  some  ^randomness”  from  run  to  tun, 
a  bias  angle  (ST)  Is  added  to  the  target  a^edt  angle 
In  SCTAMP.  It  was  suggested^py  R;  Si^lcl^^that  this  angle  be 
a  random  variable,  uniform  over^M,  that  Is  chosen  at  the 
beginning  of  each  run.  / 

The  changes  made  since  5  Hay  76  (the  version  described  In  MRI  Report  149-21) 
are  Indicated  by  a  double  $$  beginning  In  calumn  73.  Only  MAIM,  ETGEO,  and 
SCTAMP  are  affected. ^  / 
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PROGRAM  MAIN  < INPUT. OUTPUT,  TAPE5=INPUT,  TAPE6=0UTPUT) 

THIS  IS  A  SAMPLE  MAIN  PROGRAM  FOR  TEST  PURPOSES  ONLY. 

THIS  EXTENDED-TARGET  SIMULATION  PACKAGE  HAS  BEEN  PREPARED  BY 
RL  MITCHELL  OF  MARK  RESOURCES, INC  (213-822-4955),  UNDER  CONTRACT  TO 
MIRADCOM.  IT  IS  WRITTEN  WITH  THE  INTENTION  THAT  IT  WILL  BE  MADE  PART 
OF  A  REAL-TIME  SIMULATION  PROGRAM.  ALTHOUGH  SOME  OF  THE  CODE  IS  NOT 
WRITTEN  IN  COMPLETELY  OPTIMUM  FORM  (IT  IS  MORE  EASILY  UNDERSTOOD  THIS 
WAY,  AND  THE  REVISIONS  ARE  EASILY  MADE). 

ALL  ARRAYS  IN  COMMON  SHOULD  BE  DIMENSIONED  IN  THE  MAIN  PROGRAM. 


SEE  THE  SUBROUTINES  FOR  A  DEFINITION  OF  THE  VARIABLES. 
RULES  FOR  DIMENSIONING  ARRAYS . 


X,  Y,  2 . NSCAT 

AMP.  DA.  DB.  DC,  DAD _ NSCAT  (MAYBE  SMALLER) 

TWARAY . 4*NARAY 

VR.  VI.DAZ.DEL . NTAP 

THE  CHANGES  MADE  ON  5  MAY  78  ARE  INDICATED  WITH  A  *  SIGN  IN  COLUMN  73 
THE  CHANGES  MADE  SINCE  5  MAY  78  ARE  INDICATED  WITH  A  PAIR  OF  SIGNS 


COMMON  /Tl/ 
COMMON  /T2/ 
COMMON  /T3/ 
COMMON  /T4/ 
COMMON  /T5/ 
COMMON  /T6/ 
COMMON  ni/ 
COMMON  /T8/ 


XO,  YO,  20,  XOD.  YOD,  ZOD,  PSI,  THETA,  PHI,  BP,  BQ,  BR 

CPSI,  SPSI.  CTHETA,  STHETA,  CPHI,  SPHI 

XM,  YM,  ZM,  XMD,  YMD,  ZMD 

NSCAT, ST, AMPMIN,  X(20),  Y(20),  Z(20) 

NTAP, DRTAP,  DRGATE,  XL,  PTGDSQ 
NARAY, TWARAY (404) 

NS, RO, Rl. ROD. AMP(20),  DA( 20) ,  DB (20) ,  DC (20) .  DAD(20) 
PEFF,  VR(8),  VKS).  DAZ(e),  DEL(B) 


% 


DEFINE  VARIABLES 


DATA 

DATA 

LR.  LW/5,  6/ 

NTAP/8/,  DRTAP/30.  /,  DROATE/60.  /,  XL/.  02/,  AMPMIN/ 1.  E-10/. 

DATA 

DATA 

DATA 

PTGDSQ/ 1. / 

DTIME/1. / 

NARAY/ 101/ 

ST/0.  / 

READ 

(LR.  100) 

XO.  YO,  ZO 

READ 

(LR, 100) 

XOD. YOD, ZOD 

READ 

(LR, 100) 

PSI, THETA, PHI 

READ 

(LR,  100) 

BP,  BQ, BR 

% 

READ 

(LR, 100) 

XM,  YM,  ZM 

READ 

(LR, 100) 

XMD.  YMD,  ZMD 

READ 

(LR, 101) 

NSCAT 

o  o  o 


READ  <LR. 100)  (X<K). Y<K)* Z(K>» K>1, NSCAT) 

WRITE  (LW. 200)  XO. Y0> ZO 
WRITE  (LW>  201 )  XOO. YOD, ZOD 
WRITE  <LW, 202)  PSI,  THETA.  PHI 
WRITE  (LW. 203)  BP. QQ.  QR 
WRITE  (LW. 204)  XM. YM.  ZM 
WRITE  (LW. 209)  XHD. YHD.  ZMD 
WRITE  (LW. 206)  NSCAT 

WRITE  (LW.  207)  (X(K).  Y(K). Z(K).  K>1.  NSCAT) 

SUBROUTINES  TO  BE  CALLED  FROM  MAIN  OR  DRIVER  PROGRAM . 

CALL  DATAIN 
CALL  TAPSET 
CALL  ETOEO 
CALL  ETODM(DTIME) 

WRITE  (LW.  208)  RO.  Rl.ROD 

WRITE  (LW. 209)  <AMP(K).  OA(K).  DB(K).  DC(K).  DAO(K>.  K-1.  NS) 
WRITE  (LW. 210)  DTIME 

WRITE  (LW.211)  (VR(I).  VKD.DAZd).  DEL(I).  I-l.NTAP) 
WRITE  (LW.  212)  PEFF 

STOP 

100  FORMAT  (3F 10.  1) 

101  FORMAT (19) 

200  FORMAT (//14H  XO. YO. ZO . / (20X3F12.  6) ) 

201  FORMAT (//17H  XOD.  YOD. ZOD . /(20X3F12.  6) ) 

202  F0RMAT(//19H  PSI. THETA.  PHI . /(20X3F12.  6) ) 

203  F0RMAT(//14H  BP. BQ. BR . /(20X3F12.  6) ) 

204  F0RMAT(//14H  XM. YM.  ZM . / (20X3F12.  6) ) 

209  F0RMAT(//17H  XMD. YMD.  ZMD . / (20X3F12.  6) ) 

206  F0RMAT(//11H  NSCAT . /20X112) 

207  F0RMAT(//11H  X.  Y.  Z . / (20X3F12.  6) ) 

208  F0RMAT(//19H  RO. Rl.ROD . /(20X3F12.  6) ) 

209  F0RMAT(//22H  AMP.  DA.  DB.  DC.  DAD . / (20X9F12.  6) ) 

210  F0RMAT(//11H  DTIME . /20XF12.  6) 

211  FORMAT (//19H  VR. VI. DAZ.  DEL . / (20X4F12.  6) ) 

212  F0RMAT(//10H  PEFF . /20XE12.  9) 

END 
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SUBROUTINE  ETGEO 


TRANSFORMATION  TO  RADAR  SPACE  FOR  N-POINT  SCATTERER  MODEL 

IN  THIS  SUBROUTINE  WE  BEGIN  WITH  THE  MODEL  OF  AN  EXTENDED  TARGET  AND 

THE  ENGAGEMENT  GEOMETRY  IN  ORDER  TO  COMPUTE  THE  AMPLITUDE  AND  RADAR 

COORDINATES  FOR  EACH  SCATTERER  IN  THE  MODEL. 

THE  MODEL  IMPLEMENTED  IS  THE  SO-CALLED  MEDIUM-RANGE  MODEL  (SEE  MRI 

REPORT  132-44). 

ASSUMPTIONS  AND  LIMITATIONS . 

1.  ALL  SCATTERERS  ASSUMED  TO  BE  ILLUMINATED  BY  SAME  TRANSMIT 
ANTENNA  GAIN. 

2.  TARGET  ASSUMED  TO  BE  WITHIN  LINEAR  REGION  OF  MONOPULSE 
RECEIVE  BEAM. 

3.  THE  DOPPLER  SHIFT  OF  THE  TARGET  CG  IS  IMPLEMENTED  BY  MEANS 
OF  A  FINELY-CONTROLLABLE  DELAY  LINE  (THE  LASER  DEVICE), 

PLUS  THE  USE  OF  THE  FREQUENCY  SYNTHESIZER. 

4.  ONLY  ONE  PHYSICAL  TARGET  IS  SIMULATED  PER  CALL. 

ALL  COMMUNICATION  TO  AND  FROM  THIS  SUBROUTINE  IS  THRU  COMMON. 

ON  INPUT . 

/Tl/  XO,  YO,  20  »  TARGET  CO  IN  INERTIAL  COORDINATES 

XOD,  YOD.  ZOD  =  TARGET  CG  RATE  IN  INERTIAL  COORDINATES 

PSI, THETA. PHI  =  YAW,  PITCH.  ROLL  ANGLES 

BP,  BQ,  BR  «  YAW,  PITCH,  ROLL  ANGLE  BODY  RATES  ♦ 

/T3/  XM.  YM,  ZM  =  MISSILE  CG  IN  INERTIAL  COORDINATES 

XMD,  YMD.  ZMD  >  MISSILE  CG  RATE  IN  INERTIAL  COORDINATES 

/T4/  NSCAT  »  NUMBER  OF  SCATTERERS  IN  TARGET  MODEL 
ST  «  APPROXIMATE  PHYSICAL  SIZE  OF  TARGET 
AMPMIN  «  AMPLITUDE  THRESHOLD  FDR  SCATTERERS 
X,  Y,  Z  «  ARRAYS  CONTAINING  SCATTERER  LOCATIONS  IN  TARGET 
COORDINATES 

/TS/  NTAP  «  NUMBER  OF  TAPS  IN  TAPPED  DELAY  LINE 
DRTAP  a  SPACING  BETWEEN  TAPS  (2-WAY  RANGE) 

014  OUTPUT . 

/T2/  CPSI.SPSI,  .  .  .  ETC  «  SINES  AND  COSINES  OF  TARGET  ANGLES 

/T7/  NS  -  NUMBER  OF  SCATTERERS  VISIBLE 
RO  •  RANGE  TO  TARGET  CG 

R1  •  RANGE  TO  FIRST  TAP 
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ROD  -  RANCE  RATE  OF  TARGET  CG 

AnP(J}  -  AMPLITUDE  OF  J-TH  SCATTERER 

DA(J)  »  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 

DB(J)  “  INCREMENTAL  B-VECTOR  OF  J-TH  SCATTERER 

DC(J)  =  INCREMENTAL  C-VECTOR  OF  J-TH  SCATTERER 

DAO(J)  »  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

T»e  TARGET  CG  AND  MISSILE  CG  COORDINATES  ARE  IN  AN  INERTIAL  COORDINATE 
SYSTEM  REFERENCED  TO  THE  GROUND  (XY-PLANE  PARALLEL  TO  GROUND.  Z-  DOWN) 

THE  ABC-VECTORS  ARE  DEFINED  AS . 

A  -  FROM  THE  TARGET  TO  THE  MISSILE 

B  -  PARALLEL  TO  THE  GROUND.  TO  THE  LEFT  AS  VIEWED  FROM  MISSILE 
C  -  PERPENDICULAR  TO  A  AND  B  IN  RIGHT-HAND  COORDINATE  SYSTEM 

THE  TARGET  COORDINATES  ARE . 

X  -  TARGET  LONGITUDINAL  AXIS,  POSITIVE  IN  DIRECTION  OF  NOSE 
Y  -  IN  DIRECTION  OF  RIGHT  WING 
Z  -  DOWN 

THE  BODY  RATES  ARE  DEFINED  AS . 

BP  -  CW  ROTATION  RATE  ABOUT  TARGET  X-AXIS 
BQ  -  CW  ROTATION  RATE  ABOUT  TARGET  Y-AXIS 
RQ  -  CW  ROTATION  RATE  ABOUT  TARGET  Z-AXIS 

THE  DIRECTION  OF  ROTATION  IS  DEFINED  LOOKING  OUT  FROM  THE  COORDINATE 
ORIGIN. 

SEE  SUBROUTINE  XFORM  FOR  A  DEFINITION  OF  THE  YAW.  PITCH.  AND  ROLL 
ANGLES. 

THE  RFSS  CHAMBER  COORDINATES  ARE  ASSUMED  TO  BE  PARALLEL  TO  THE  ABC- 
VECTORS.  RANGE  IS  IN  -A  DIRECTION.  RIGHT  AZIMUTH  IN  -B  DIRECTION.  AND 
UP  ELEVATION  IN  -C  DIRECTION. 

ALL  DISTANCES  (INCLUDING  WAVELENGTH)  MUST  BE  IN  THE  SAME  UNITS.  ALL 
ANGLES  MUST  BE  IN  RADIANS. 

DIMENSION  A(3).  B<3).C(3). WA(3).  WB(3).  WC(3) 

COMMON  /Tl/  XO.  YO.ZO.XOD.YOD.  ZOD.PSI.  THETA.  PHI.  BP.  BQ.BR  « 

COMMON  /T2/  CPSI.  SPSI.  CTHETA,  STHETA.  CPHI.  SPHI 

COMMON  /T3/  XM.  YM.  ZM.  XMD,  YMD.  ZMD 

COMMON  /T4/  NSCAT,  ST,  AMPMIN.  X(20),  Y(20),Z<20) 

COMMON  /T5/  NTAP, DRTAP 

COMMON  mi  NS.  RO.  Ri.ROD.  AMP(20>.DA(20).0B(20},DC(20).DA0(20) 
COMPUTE  SINES  AND  COSINES  OF  ANGLES 


ooo  nnnoono  ooo  non  oon  nnn  ooo 


C 


CALL  SINCOSCPSI. SPSl.  CPSZ) 

CALL  S I NCOS < THETA.  STHETA,  CTHETA) 

CALL  SINCOS(PHZ. SPHI.  CPHI) 

COMPUTE  RANGE  TO  TARGET  CG  AND  A-VECTOR 

A(l)«XM-XO 

A(2>*YM-Y0 

A(3)-ZM-Z0 

RO»SQRT  ( A  ( 1 )  **2-»-A  ( 2 )  «»2-<<A  ( 3 )  *«2  > 

A(1)-A(1)/R0 

A(2)-A(2>/R0 

A(3)«A(3)/R0 

COMPUTE  RANGE  TO  FZRST  TAP  ^ 

R1-RO-.  5*(NTAP-1)*DRTAP 

COMPUTE  RANGE  RATE  OF  TARGET  CG 

ROO«A  < 1 ) ♦<  XOD-XMD) +A<2) ♦ <  YOD-YMD) +A ( 3) * ( ZOD-ZMD ) 

COMPUTE  B-  AND  C-VECTORS 

RHO»SaRT  <  A  (  1 )  »«2'»-A  ( 2 )  *«2 ) 

B(1)«-A(2)/RH0 
B(2)>  A(l)/RHO 
B(3)«0. 

C(1)»A<3)«B<2) 

C(2)>  A<3)«B<1) 

C(3)-RH0 

TRANSFORM  A-,  B-,  AND  C-VECTORS  TO  TARGET  COORDZNATES 

CALL  XFORM<A.  WA) 

CALL  XFORM(B.WB) 

CALL  XFORM(C.WC) 


COMPUTE  TARGET  ASPECT  ANGLE  (ALPHA-AZZMUTH.  BETA-ELEVATZON.  tt 

ANOL-ANGLE  TO  ROLL  AXZS) 

ALPHA>ATAN2(WA(2).WA(1))  «« 

8BETA— WA(3)  %% 

BETA«ATAN2 ( SBETA. SORT ( 1 .  -SBETA*«2 ) ) 

AN0L«ATAN2(SQRT(1. -WA(1)«*2).WA(1)}  «« 


LOOP  OVER  SCATTERERS 
L«1 


DO  20  Ksl.NSCAT 
SAMP^SCTAMP ( K.  AN6L ) 

ZFCSAMP.  LE.  AMPMIN)  GO  TO  20 
AMP(L)>SAnP 
C 

C  COMPUTE  INCREMENTAL  A> B. C  COORDINATE 
C 

DA(L>-X(K)«WA<  1  )-»-Y(K)«MA(2)-4-Z(K)«WA<3) 
DB(L)-X(K)«WB<l)-«-Y<K)«NB(2)-i-Z(K>«UD(3) 
DC(L)aX(K)*UC(l)+Y(K)«WC(2)-<-Z(K>*WC(3} 

C 

C  COMPUTE  INCREMENTAL  A-VECTOR  RATE  (SMALL  ANGLES  ARE  ASSUMED) 
C 

XKD-  Z(K)«BQ-Y(K)«BR 
YKD—Z  ( K  )  *BP+X  <  K )  *BR 
ZKD-  Y(K)«BP-X<K)«BQ 
DAD  (  L )  «XKD«WA  <  1 )  -»-YKD«MA  (  2  )  '•■ZKD«WA  (  3 ) 

L-L+1 

20  CONTINUE 
N8-L-1 
RETURN 
END 


ooonnnonnonnnnoonnnoonnnnononnoonoononoono 


SUBROUTINE  ETGDM<DTinE) 

OLXNT  AND  DOPPLER  MODULATION  FOR  N-POINT  SCATTER  MODEL 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  GLINT  OFFSETS  AND  MODULATION  SIGNALS 
APPLIED  TO  EACK  TAP  OF  THE  TAPPED-DELAY  LINE.  IT  IS  TO  BE  CALLED 
AFTER  ETGEO  TRANSFORMS  COORDINATES  TO  RADAR  SPACE.  IT  WILL  USUALLY 
BE  CALLED  MORE  FREQUENTLY  THAN  ETGEO.  IT  IS  ALSO  THE  BEST  SUBROUTINE 
TO  PLACE  IN  THE  AP120B. 

EXCEPT  FOR  TIME,  ALL  COMMUNICATION  TO  AND  FROM  THIS  SUBROUTINE  IS  THRU 
COMMON. 

ON  INPUT . 

DTIME  «  TIME  SINCE  LAST  UPDATE  IN  TARGEO 

/T5/  NTAP  «  NUMBER  OF  TAPS  IN  TAPPED  DELAY  LINE 
DRTAP  «  SPACING  BETWEEN  TAPS  (2-WAY  RANGE) 

XL  *  WAVELENGTH 

PTGDSQ  «  PRODUCT  OF  TRANSMIT  POWER,  GAIN,  AND  SQUARE  OF 
RFSS  CHAMBER  LENGTH 

/T7/  NS  -  NUMBER  OF  SCATTERERS  VISIBLE 

RO  »  RANGE  TO  TARGET  CG 

R1  *  RANGE  TO  FIRST  TAP 

AMP<J)  ■  AMPLITUDE  OF  J-TH  SCATTERER 
DA(J)  »  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 
DB(J)  *  INCREMENTAL  B-VECTOR  OF  J-TH  SCATTERER 
DAO(J)  »  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

ON  OUTPUT . 

/TS/  PEFF  »  EFFECTIVE  RADIATED  POWER  AT  RFSS  ARRAY 

VR<I)  e  IN-PHASE  MODULATION  SIGNAL  TO  I-TH  TAP 
VKI)  a  QUADRATURE  MODULATION  SIGNAL  TO  I-TH  TAP 
DAZ(I)  ■  GLINT  OFFSET  (AZIMUTH)  FOR  I-TH  TAP 
DEL(I)  •  GLINT  OFFSET  (ELEVATION)  FOR  I-TH  TAP 

THE  PARAMETER  PMIN  IS  JUST  SOME  SMALL  NUMBER  TO  PREVENT  DIVIDE  BY  ZERO 

ARRAYS  VBR, VBI,  VCR, VCI  MUST  BE  DIMENSIONED  AS  LARGE  AS  NTAP 

DIMENSION  VBR(8),  VBKS),  VCR(S),  VCKS) 

DIMENSION  SS(4),CC(4) 

DIMENSION  TW(4) 

COMMON  /T5/  NTAP,  DRTAP, DRG ATE,  XL,  PTGDSQ 

COMMON  /T7/  NS,  RO,  R1,R0D,AMP(20),DA(20),DB(20),DC(20),DAD(20) 
COMMON  /TS/  PEFF,  VR(S),  VI(8),DAZ(8),DEL(8) 

DATA  PMIN/1,  E-IO/ 


DATA  FOURPI/12.  5663706/ 

C 

C  ZERO  ARRAYS  ^ 

C 

CALL  XMIT<-NTAP,0.  .  VR) 

CALL  XMIT(-NTAP,  0.  .  VI ) 

CALL  XMIT(-NTAP* 0.  .  VBR) 

CALL  XMIT(-NTAP,  0.  .  VBI ) 

CALL  XMIT < -NTAP,  0.  >  VCR ) 

CALL  X«IT<-NTAP,  0.  .  VCI ) 

CALL  XMIT<-NTAP,0.  ,DAZ) 

CALL  XM1T(-NTAP,0.  ,DEL) 

C 

C  LOOP  OVER  NS  SCATTERERS 
C 

DO  40  Jal.NS 
C 

C  COMPUTE  TAP  NUMBER  OF  FIRST  TAP  (ITAP)  AND  FRACTION  (P) 
C 

R-RO- <  DA  ( J )  4-DAD  ( J )  «DT  I  ME  ) 

P- <  R-R 1 ) /DRTAP4- 1 00. 

XTAP-P 

P-P-ITAP 

XTAP»ITAP-100 

C 

C  COMPUTE  RANGE  DIFFERENCE  FROM  TAP  NUMBER  ITAP 
C 

DR-<P4-1.  )*DRTAP 
C 

C  FIND  TAP  WEIGHTS 
C 

CALL  TAPWTS<P»TW) 

C 

C  COMPUTE  PHASE  ON  FOUR  TAPS 
C 

DO  20  X-1,4 

CALL  SXNCOS(-FOURPI*DR/XL»S»C) 

88( X )-S«AMP ( J } «TW ( I ) 

CC ( X } «C«AMP ( J ) *TW  < I ) 

DR>DR-DRTAP 
20  CONTXNUE 
C 

C  LOOP  OVER  UP  TO  FOUR  TAPS  AND  INCREMENT  ARRAYS 
C 

XFdTAP.  GT.  NTAP)  GO  TO  40 
XFdTAP.  LT. -2)  GO  TO  40 
Xi-MAXO<ITAP. 1) 

X2-MXN0<  XTAP4-3,  NTAP) 

XX-Xl-XTAP 
DO  30  X-Ii< 12 


IX-XI+1 

VR  (I)»VR  (I)+CC(H) 

VX  (X>-VI  (I)+SS<H) 

VBR(X>-VBR(I)-»-CC(II)*DB(J) 

VBX(I)«VBI(I>+SS<II)«DB(J) 

VCR<I)«VCR<I)+CC<H)*DC<J) 

VCX(l)»VCI(I}-t-SS(Il)«OC<J) 

30  CONTINUE 
40  CONTINUE 
C 

C  COMPUTE  OLINT  OFFSETS  FOR  EACH  TAP  AND  PEAK  POWER 
C 

PEAK-O. 

DO  90  Xsl.NTAP 
POU*VR ( I ) **2+VI ( I ) *«2 
XF(POU.  GT.  PEAK)  PEAK-POW 
XF<POW.  LT.  PMIN)  GO  TO  50 

OAZ(I)  — <VDR<I)*VR<I)-t'VBI(I)«VI(I>>/(RO«POW) 
DEL(I)  — (VCR(I)*VR<I)+VCI(I>*VI(I>)/(RO*POW> 
90  CONTINUE 
C 

C  NORMALIZE  AMPLITUDE 
C 

ANORM«SQRT(PEAK) 

DO  60  I*liNTAP 
VR(X)«VR<I)/ANORM 
VX(I)«VI<I)/ANORM 
60  CONTINUE 
C 

C  COMPUTE  EFFECTIVE  RF  POWER 
C 

PEFF«PEAK*PTGDSQ/ ( FOURP I »R0»»4  > 

RETURN 

END 
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SUBROUTINE  XFORM<A.  U) 

IN  THIS  SUBROUTINE  WE  TRANSFORM  A  VECTOR  <A)  IN  INERTIAL  COORDINATES 
TO  A  VECTOR  <W>  IN  TARGET  COORDINATES.  THE  COORDINATE  ROTATIONS,  IN 
THE  ORDER  OF  APPLICATION.  ARE . 

PSI  -  CW  ROTATION  OF  Z-AXIS  (YAW) 

THETA  «  CW  ROTATION  OF  Y-AXIS  (PITCH) 

PHI  «  CW  ROTATION  OF  X-AXIS  (ROLL) 

THE  DIRECTION  OF  ROTATION  IS  DEFINED  LOOKING  INTO  THE  COORDINATE 
ORIGIN.  IN  THIS  SUBROUTINE  THE  SINES  AND  COSINES  OF  THE  ANGLES  ARE 
INPUT  THRU  COMMON  /T2/. 

DIMENSION  A(3),W(3) 

COMMON  /T2/  CPSI.SPSI.CTHETA.  STHETA,  CPHI.SPHI 
UX«A ( 1 ) *CPSI -A ( 2 ) *SPS I 
UY*A(1)«SPSI-^A(2)*CPSI 
UZ-A(3) 

VX-  UX*CTHETA+UZ*STHETA 
VY«  UY 

VZ«-UX*STHETA+UZ*CTHETA 

W(1)«VX 

W<2)-VY*CPHI-VZ*SPHI 

W<3)»VY*SPHI+VZ*CPHI 

RETURN 

END 
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SUBROUTINE  TAPUTS<P.  TU) 

IN  THIS  SUBROUTINE  FOUR  TAP  WEIGHTS  ARE  RETURNED  IN  ARRAY  TU  ACCORDING 
TO  THE  FRACTION  P.  THE  WEIGHTS  ARE  EXTRACTED  FROM  A  PRECOMPUTED  TABLE 
(SEE  SUBROUTINE  TAPSET). 

ARRAY  TWARAY  IS  USED  AS  IF  IT  WERE  DIMENSIONED  (4.  NARAY). 

DIMENSION  TW(4) 

COMMON  nU  NARAY.  TWARAY(l) 

DATA  LW/6/ 

1NDEX>(NARAY-1)*P+1.  9 

CALL  XMIT(4.  TWARAY(4«IN0EX-3).  TW) 

>6.  2iF5X02/ . WT.P  HOKTAMROF  C 

UT.P  >002.  WL(  ETIRW 

RETURN 

END 
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SUBROUTINE  TAPSET 

IN  THIS  SUBROUTINE  THE  TAP  WEIGHT  TABLE  IS  COMPUTED.  IT  IS  A 
COMPANION  SUBROUTINE  TO  TAPWTS.  AND  IT  IS  TO  BE  CALLED  AS  AN  INITIAL¬ 
IZATION  STEP  PRIOR  TO  THE  BEGINNING  OF  THE  SIMULATED  MISSION. 

/T5/  DRTAP  >  SPACING  BETWEEN  TAPS  (2-WAY  RANGE) 

DRGATE  >  SPACING  BETWEEN  RECEIVER  GATES  (2-WAY  RANGE) 

ARRAY  TWARAY  MUST  BE  DIMENSIONED  AS  LARGE  AS  4«NARAY. 

DIMENSION  A(4.4)<X(4) 

COMMON  /T5/  NTAP.  DRTAP» DRGATE 
COMMON  nU  NARAY.  TWARAY(l) 

D-DRTAP/DRGATE 

L-1 

DO  30  K»l.  NARAY 
P-(K-1 ) /FLOAT (NARAY-1 ) 

DO  10  J»l»4 
X(J)-CHI(D*(P+2-J)) 

10  CONTINUE 
DO  20  I-1.4 
DO  20  J>1»4 
A(I.  J)-CHI(D«(I-J) ) 

20  CONTINUE 

CALL  SIMa(A<  X>  4# lERR) 

IFdERR.  GT.  0)  STOP 
CALL  XMIT ( 4.  X«  TWARAY (D) 

L*L+4 

30  CONTINUE 
RETURN 
END 
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FUNCTION  CHKP) 

RANGE  GATE  RESPONSE.  THE  ARGUMENT  P  IS  THE  RANGE  MISMATCH  NORMALIZED 
TO  THE  RECEIVER  GATE  SPACING.  INTERPOLATION  IS  USED  ON  THE  SAMPLES 
STORED  IN  THE  A-ARRAY.  WHERE  THE  SPACING  IS  0.  1  UNIT. 

THE  RESIDUAL  ERROR  IN  THE  INTERPOLATION  IS  LESS  THAN  .  0003 

P  MUST  BE  LESS  THAN  1.  5  IN  MAGNITUDE. 

THE  SAMPLES  ARE  OF  THE  RESPONSE  DERIVED  IN  MRI  REPORT  149-4. 

DIMENSION  A(18) 

DATA  A/ 1.  00000.  .98104.  .92193.  .81903.  .67431.  .50112.  .32385. 

1  .17071.  .06308.  .  00731. -.  00651.  .00182.  .01262.  .01458. 

2  .  00713. -.  00313. -.  00898. -.  00762  / 

H*10.  •ABS(P) 

1F<H.  GT.  15.  )  STOP  55 

I-H 

H^-I 

IPl-1+1 

1P2-1+2 

IP3-I+3 

IFCI.LE.  0)  1»2 

CHI*-.  166667*H*<H-1.  )*<H-2.  )*A(I>+.  5*<H**2-1.  )*<H-2.  )*A(IP1> 

1  -.  5*H*<H-»-l.  )*<H-2.  )*A(IP2)+.  166667*H*(H**2-1.  )*A<IP3) 

RETURN 

END 
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SUBROUTINE  SinQ(A.B.N.ZERR) 

SOLVES  SET  OF  N  SIMULTANEOUS  EQUATIONS . 

A  *  X  -  B  SUM  (A(l.  J>*X(J))  -  B(I) 

WHERE  ARRAY  A  IS  2-DIMENSIONAL.  ARRAY  X  IS  RETURNED  IN  ARRAY  B,  AND 
ARRAY  A  IS  DESTROYED.  COMPUTATION  IS  VALID  IF  lERR^O 

DIMENSION  A(1).B<1) 
lERR  -  0 

IF  (N.  6T.  0)  GO  TO  10 

ZERR  *  1 

RETURN 

FORWARD  SOLUTION 

10  TOL  -  0.  0 
K8  -  0 
JJ  «  -N 
DO  69  J  -  1«N 
UY  -  J  +  1 
JJ  -  JJ  4-  N  1 
BIOA  -  0. 

IT  -  JJ  -  J 
DO  30  I  «  J«N 

SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 
IJ  -  IT  +  I 

IF  (ABS(BIGA)-ABS(A(ZJ)))  20,  30,  30 
20  BIOA  s  A(IJ) 

ZMAX  »  I 
30  CONTINUE 

TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

IF  (ABS(BIOA)-TOL)  39*  39.  40 

39  ZERR  >  2 
RETURN 

INTERCHANGE  ROWS  IF  NECESSARY 

40  11  -  J  N*(J-2) 

IT  -  IMAX  -  J 
DO  90  K  -  J.N 

XI  -  11  N 
X2  «  11  XT 
SAVE  -  A(X1) 

A(X1)  -  A<X2) 


f»oo  non  noo 


A(X2)  -  SAVE 

DIVIDE  EQUATION  BY  LEADING  COEFFICIENT 

90  A(I1)  -  A<I1)/BIGA 
SAVE  -  B(IMAX) 

BdNAX)  -  B(J) 

B<J)  -  SAVE/BIGA 

ELIMINATE  NEXT  VARIABLE 

IF  (J-N)  95.  70.  99 

99  IQS  -  N«(J-1) 

DO  69  IX  -  JY.N 
IXJ  «  IQS  *  IX 
IT  ■  J  -  IX 
DO  60  JX  «  JY.N 
IXJX  -  N«(JX-1)  +  IX 
JJX  -  IXJX  +  IT 

60  AdXJX)  -  A(IXJX)  >  (A<IXJ)«A(JJX)) 

69  B(IX)  «  B(IX)  -  (B(J)«A(IXJ)) 

BA  SOLUTION 

70  NY  •  N  -  1 
IT  •  N«N 

DO  60  J  «  l.NY 
lA  -  IT  -  J 
IB  »  N  -  J 
IC  «  N 

DO  80  K  -  l.J 

B(IB)  «  B<IB)  -  A(IA)«B(IC) 

XA  ■  XA  -  N 
80  IC  -  IC  -  1 
RETURN 
END 
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SUBROUTINE  XmT(N.  A.  B) 

IN  THIS  SUBROUTINE  WE  EITHER  TRANSMIT  ARRAY  A  TO  ARRAY  B  (IF  N.  OT.  0) 
OR  WE  TRANSMIT  THE  CONSTANT  A  TO  ARRAY  B  (IF  N.  LT.  0).  IN  EITHER  CASE 
THE  ARRAY  LENGTH  IS  lABS(N). 

THIS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEMBLY  LANGUAGE 

DIMENSION  A(1),B(1) 

1F(N)  10i20.  25 
10  NN— N 
"  AA-A(l) 

DO  15  K-l.NN 
B(K)-AA 
IS  CONTINUE 
20  RETURN 
25  DO  30  K-l.N 
B(K)«A(K) 

90  CONTINUE 
RETURN 
END 
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SUBROUTINE  SINCOS(ARe<  S,  C) 

THXS  SUBROUTINE  SHOULD  BE  WRITTEN  IN  ASSEMBLY  LANGUAGE*  USING  THE 
TABLE-LOOKUP  METHOD  DESCRIBED  BY  MITCHELL  (RADAR  SIGNAL  SIMULATION). 


S-SIN(ARG) 

C«COS(ARG} 

RETURN 

END 
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SUBROUTINE  ETGDKDTIME} 


GLINT  AND  DOPPLER  MODULATION  FOR  N-POINT  SCATTER  MODEL 
NO  RANGE  EXTENSION 
SUBROUTINE  REPLACES  ETGDM 
ON  INPUT . 

DTIME  -  TIME  SINCE  LAST  UPDATE  IN  TARGEO 
/T5/  XL  -  WAVELENGTH 

PTGDSQ  -  PRODUCT  OF  TRANSMIT  POWER.  GAIN.  AND  SQUARE  OF 
RFSS  CHAMBER  LENGTH 

/T7/  NS  -  NUMBER  OF  8CATTERERS  VISIBLE 

AMP(J)  -  AMPLITUDE  OF  J-TH  SCATTERER 
DA<J)  -  INCREMENTAL  A-VECTOR  OF  J-TH  SCATTERER 

DB<J)  -  INCREMENTAL  B-VECTOR  OF  J-TH  SCATTERER 

DC<J)  •  INCREMENTAL  C-VECTOR  OF  J-TH  SCATTERER 

DAD<J>  -  INCREMENTAL  A-VECTOR  RATE  OF  J-TH  SCATTERER 

ON  OUTPUT . 

/TO/  PEFF  -  EFFECTIVE  RADIATED  POWER  AT  RFSS  ARRAY 

/T9/  VR.VI  «  DOPPLER  MODULATION  SIGNAL 

DR.DAZ.  DEL  -  RANGE.  AZIMUTH.  AND  ELEVATION  GLINT  OFFSETS 

COMMON  /T5/  NTAP. DRTAP.  DRGATE.  XL.  PTGDSQ 

COMMON  /T7/  NS.  RO.  R1.R0D.AMP(20).DA(20).DB(20>.DC(20>.DAD(20> 

COMMON  /T8/  PEFF 

COMMON  /T9/  VR.VI.  DR.  DAZ.  DEL 

DATA  FOURPI/12.  5663706/ 

ZERO  ACCUMULATORS 


VR«0. 

V1«0. 

VAR-0. 

VAI-0. 

VBR-0. 

VBI-O. 

VCR-0. 

VCI-0. 

LOOP  OVER  NS  8CATTERER8 


DO  40  J-1.N8 
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CALL  8INC0S(F0URPI*(DA(J)-*-DAD(J)*DTIME)/XL*8.C> 
C-C«AMP(J) 

8>8«AMP<J) 

VR  -VR  -^C 
VI  -VI  +8 
var-var-k:*da(J) 

VAI«VA1-^S«DA(J> 

VBR-VBR-<'C«DB(J) 

VBI-VBI'i-S«OB(J) 

VCR-VCR-«-C«OC(J) 

VCI-VCI-t‘8«0C(J) 

40  CONTINUE 

P0U-VR*«2**>VI*«2 

ANPL-8aRT(P0U) 

CONPUTE  GLINT  0FF8ET8 

OR  ->(VAR«VR-*-VAI«VI)/POW 
DAZ— (  VBR«VR4-VBI«VI  )  /  (RO«PaU) 

DEL— <  VCR*VR+VC  I *VI )  /  <  RO*POW ) 

COMPUTE  EFFECTIVE  RF  POWER 

PEFF-P0U*PT0D8a/ <  FOURP I •R0««4 ) 

NORMALIZE 

VR-VR/AMPL 

VI-VI/AMPL 

RETURN 

END  . 
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SUBROUTINE  DATAIN 

READS  TARGET  SCATTERING  DATA  SUPPLIED  BY  M.  MUMFORD  (SEE  SCTAMP). 

DIMENSION  IA(1}< AA(4)i  XX(4>.yY(4)>  ZZ(4) 

COMMON  /DP/  P(IOO). IP(IOO) 

COMMON  /DQ/  Q(918) 

COMMON  /T4/  NSCAT 
DATA  LR.  LU/5i6/ 

NSCAT-10 

DO  20  I»l.  NSCAT 
PRINT  99.  I 
L-10*(l-l) 
to  L-L-t-l 

READ  (LR.  100)  lA.  P(L).  AA.  XX.  YY.  ZZ 
WRITE  (LW.  100)  lA.  P(L).  AA.  XX.  YY.  ZZ 
IP(L)«M 
ZA«IA-2 

CALL  XMIT(17.  lA.  Q(M)) 

M-N+17 

1P(P(L).  LT.  180.  )  GO  TO  10 
SO  CONTINUE 
RETURN 

99  FORMAT (/29H  TARGET  DATA  FOR  SCATTERER  N0I3//) 

100  PORMATdXIl.  12XF8.  3.  4E14.  8/ (22X4E14.  8)) 

END 
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FUNCTION  SCTAMP(K.  AN6L)  ** 

IN  THIS  SUBROUTINE  WE  COMPUTE  THE  AMPLITUDE  <SQRT(RCS))  OF  THE  K-TH 
SCATTERER  AS  VIEWED  FROM  THE  TARGET  ASPECT . 

ANOL  B  ANGLE  FROM  ROLL  AXIS  MEASURED  FROM  NOSE  <RAD)  ** 

IN  ADDITION  IN  COMMON  /T4/ .  M 

ST  -  BIAS  THAT  IS  ADDED  TO  ANGL  <RAD)  «« 

«« 


THIS  SUBROUTINE  ACCESSES  TARGET  DATA  SUPPLIED  BY  MIKE  MUMFORD  AT  NWC/ 
CHINA  LAKE  IN  THE  FORMAT  DEFINED  BY  A  COMPUTER  PROGRAM  WRITTEN  5/11/78 
BY  E.  HUTTON  X3219. 


DIMENSION  IA(1). AA(4). XX(4),  YY(4),  ZZ(4} 

COMMON  /T4/  NSCAT.  ST.AMPMIN,  X(20).  Y(20>«Z(20) 
COMMON  /DP/  P(IOO).  IPdOO) 

COMMON  /DQ/  Q(918) 


ANG-ABS  <  ANGL+ST  >  *57.  2957795 
1P(ANG.  GT.  180.  )  ANG>180. 

Il-10*<K-i)+l 
I2»Il-*-9 
DO  20  I-Il.  12 
1F(ANG.  LT.  P(I) )  60  TO  25 
20  CONTINUE 
25  M»IP(I) 

CALL  XMIT(17.  Q(M)» lA) 

IFdA)  30,31.32 

30  8CTAMP-AA(1)^ANG«<AA<2)*ANG*(AA(3>'^ANG*AA(4>  >  ) 

GO  TO  35 

31  8CTAMP>EXP(AA(1)*ANG*AA<2)) 

GO  TO  35 

32  SCTAMP-EXP  <  AA  ( 1 )  •^AN6*AA  (  2  )  )  -EXP  (  AA  ( 3  >  •••ANG«AA  ( 4  )  ) 
35  SCTAMP-.  09004*SCTAMP 

-  IFCSCTAMP.  LT.  AMPMIN)  RETURN 

X(K>>XX(l>*ANG*<XX(2>*ANG*CXX(3>i-AN0*XX(4>)> 
Y(K)bYY<1>+AN0*<YY<2)+ANG*(YY(3)*ANG*YY(4) > ) 
7(K)-ZZ(1)*AN0*(ZZ(2)*AN0*(ZZ(3)*ANG*ZZ(4>)) 
Y<K)— Y<K) 

Z(K>— Z(K} 

R^TUBN  .  .  V-  k  - - 

END  -  - 


«« 


